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Abstract 



Dilepton production in pp and Au + Au nucleus-nucleus collisions at y/s = 200 GeV 
as well as in In+In and Pb+Au at 158 A-GeV is studied within the microscopic HSD 
transport approach. A comparison to the data from the PHENIX Collaboration at 
RHIC shows that standard in-medium effects of the p, to vector mesons - compatible 
with the NA60 data for In + In at 158 A-GeV and the CERES data for Pb + Au at 
158 A-GeV - do not explain the large enhancement observed in the invariant mass 
regime from 0.2 to 0.5 GeV in Au + Au collisions at yfs = 200 GeV relative to pp 
collisions. 
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While the properties of hadrons are rather well known in free space (embedded in 
a nonperturbative QCD vacuum) the masses and lifetimes of hadrons in a baryonic 
and/or mesonic environment are subject of current research in order to achieve a 
better understanding of the strong interaction and the nature of confinement. In this 
context the modification of hadron properties in nuclear matter are of fundamental 
interest (cf. Refs. pl2|3|H|f5] ) since QCD sum rules [2.3.6] as well as QCD inspired 
effective Lagrangian models [TPfTf9yT0] predict significant changes e.g. of the vector 
mesons (p, u and 0) with the nuclear density pn and/or temperature T [TTJT2"] . 




A modification of vector mesons has been seen experimentally in the enhanced pro- 
duction of lepton pairs above known sources in nucleus-nucleus collisions at Super- 
Proton-Synchroton (SPS) energies [TBTTi] . As proposed by Li, Ko, and Brown [15] 
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and Ko et al. [16] the observed enhancement in the invariant mass range 0.3 < M < 
0.7 GeV might be due to a shift of the p-meson mass following Brown/Rho scaling [1] 
or the Hatsuda and Lee sum rule prediction J2]. The microscopic transport studies 
in Refs. [TTfl7lll8|[T5] for these systems have given support for this interpretation. 
On the other hand also more conventional approaches that describe a melting of 
the p-meson in the medium due to the strong hadronic coupling (along the lines of 
Refs. pfTO] ) have been found to be compatible with the early CERES data [T2"f 171120] . 
This ambiguous situation has been clarified to some extent in 2006 by the NA60 Col- 
laboration since the invariant mass spectra for p + p~ pairs from In+In collisions at 
158 A-GeV favored the 'melting p' scenario [2T]. Also the more recent data from the 
CERES Collaboration (with enhanced mass resolution) |22j show a preference for 
the 'melting p' picture. 



In 2007 the PHENIX Collaboration has presented first dilepton data from pp and 
Au + Au collisions at Relativistic-Heavy- Ion-Collider (RHIC) energies of \fs = 200 
GeV [23] which show an even larger enhancement in Au + Au reactions (relative to 
pp collisions) in the invariant mass regime from 0.15 to 0.6 GeV than the data at 
SPS energies [2"T||2"2"] . The question arises if this sizeable enhancement might also be 
attributed to in-medium modifications of the p and u mesons as at SPS energies 
[TTp2] or if new radiative channels from the strong Quark-Gluon Plasma (sQGP) 



have been seen. 



The answer to this question is nontrivial due to the nonequilibrium nature of the 
heavy-ion reactions and covariant transport models have to be incorporated to dis- 
entangle the various sources that contribute to the final dilepton spectra seen exper- 
imentally. In this study we aim at contributing to this task employing an up-to-date 
relativistic transport model (HSD) that incorporates the relevant off-shell dynamics 
of the vector mesons. The HSD transport model [TTfl8f24j has been used for the 
description of pA and AA collisions from SIS to RHIC energies and lead to a fair re- 
production of hadron abundancies, rapidity distributions and transverse momentum 
spectra. We recall that in the HSD approach nucleons, A's, N*(1440), N*(1535), A, S 
and S* hyperons, E's, S*'s and JTs as well as their antiparticles are included on the 
baryonic side whereas the 0~ and 1~ octet states are incorporated in the mesonic 
sector. Inelastic baryon-baryon (and meson-baryon) collisions with energies above 
\fs th ~ 2.6 GeV (and ^/!hh — 2.3 GeV) are described by the Fritiof string model 
[25] whereas low energy hadron-hadron collisions are modeled in line with experi- 
mental cross sections. Low energy cross sections such as threshold meson production 
in proton-neutron (pn) collisions - which are scarcely available from experiments - are 
fixed by proton-proton (pp) cross sections and isospin factors emerging from pion- 
exchange diagrams. Since we address ultrarelativistic collisions at SPS and RHIC 
energies such 'low energy uncertainties' are of minor relevance here. As pre-hadronic 
degrees of freedom HSD includes 'effective' quarks (antiquarks) and diquarks (antidi- 
quarks) which interact with cross sections in accordance with the constituent quark 
model (cf. Refs. |26j). 
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Compared to our earlier studies in Refs. [TTf20 j a couple of extensions have been 
implemented such as 

• off-shell dynamics for vector mesons - according to Refs. [27J - and an extended 
set of vector meson spectral functions [28] 

• extension of the LUND string model to include 'modified' spectral functions for 
the hadron resonances in the string decays. 

As demonstrated in Ref. [2B] the off-shell dynamics is particularly important for 
resonances with a rather long lifetime in vacuum but strongly decreasing lifetime in 
the nuclear medium (especially u and mesons) but also proves vital for the correct 
description of dilepton decays of p mesons with masses close to the two pion decay 
threshold 0- For a detailed description of the various hadronic channels included for 
dilepton production as well as the off-shell dynamics we refer the reader to Ref. [2"8"] 
where we have focused on e + e~ production in the 1 to 2 A-GeV energy range. 

Before we step to a comparison with experimental data we show in Fig. [1] the time 
evolution of the mass distribution of p mesons for central (b=0.5 fm) In + In colli- 
sions at 160A GeV (left part) and for Au + Au collisions at 21. 3A TeV (right part) 
for the free p spectral function (upper plots) and the 'collisional broadening' scenario 
(lower plots). In the free case there are no mass components below 2m 7T whereas in 
the 'collisional broadening' scenario the p mass distribution extends down to twice 
the electron mass. The p mass distributions for times t > 15 fm/c essentially show 
the width due to the p — ► ttti decay in vacuum which is delayed at the RHIC en- 
ergy due to significantly larger Lorentz 7 factors (with respect to the calculational 
frame). In the collisional broadening case (lower plots) we observe a substantially 
larger width in the initial p meson mass distribution at the RHIC energy due to a 
larger baryon + antibaryon density for the heavier Au+Au system and a higher p 
+ meson scattering rate. The corresponding results on dilepton spectra - within the 
experimental accceptance and mass resolution - will be shown below. 

We mention that also finite temperature effects lead to a sizable broadening of the 
vector mesons spectral functions (also at baryon chemical potential ps — 0). This 
is essentially due to vector meson scattering with mesons which may contribute to 
the total width by 70-80 MeV at a temperature of ~ 170 MeV according to the 
early work by Haglin [29]. For the present study we use a simplified modeling of the 
collisional broadening by meson scattering which discards an explicit consideration 
of such 'temperature effects' in the parametrization of the vector-meson spectral 
functions (cf. Ref. [28J). However, the 'temperature effects' are partly accounted 
here by explicit meson-meson interactions which lead to a dynamical broadening in 
the vector meson mass distribution. 

Let's directly go over to the actual results from HSD for In + In collisions at 160 



This will be of particular importance in comparison to the low mass dileptons from 
PHENIX. 
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Fig. 1. Time evolution of the mass distribution of p mesons for central (b=0.5 fm) In + In 
collisions at 160A GeV (left part) and for Au + Au collisions at 21. 3A TeV (right part) for 
the free spectral function (upper plots) and for the 'collisional broadening' scenario (lower 
plots). 



A -GeV which are displayed in Fig. [2] for peripheral, semi-peripheral, semi-central 
and central collisions following the centrality definition of the NA60 Collaboration 
[2"T|30f31j : all measured events have been separated into 4 centrality bins according to 
their charged particle multiplicity dN ch /drj measured in the pseudorapidity interval 
3 < rj < 4.2 (cf. Ref. [30] for details): bin 1 (central) - 170 < dN ch /dr] < 240; bin 2 
(semi-central) - 110 < dN ch /drj < 170; bin 3 (semi-peripheral) - 30 < dN ch /dr] < 110; 
bin 4 (peripheral) - 4 < dN c h/dr) < 30. By computing dN c h/dr) within HSD for the 
same pseudorapidity interval we obtain a direct correspondence of the centrality 
bins with the impact parameter intervals: bin 1 (central) - b < 3.5 fm; bin 2 (semi- 
central) - 3.5 < b < 5.5 fm; bin 3 (semi-peripheral) - 5.5 < b < 8.5 fm; bin 4 
(peripheral) - b > 8.5 fm which we used as centrality criteria for our calculations. 
We note, furthermore, that the experimental data in Fig. [2] correspond to the 'excess 
mass spectra', which are extracted from the measured dilepton yields by subtracting 
the 'cocktail' contribution from the r],u,cf) decays [2Tj . Such a procedure allows to 
'separate' the contribution of the p mesons which is dominant in the mass region 
0.25 < M < 0.9 GeVc 2 over other dilepton sources. Following Ref. [21] we compare 
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Fig. 2. The HSD results for the mass differential dilepton spectra from the direct p meson 
decay in case of In + In at 158 A-GeV for peripheral, semi-peripheral, semi-central and 
central collisions in comparison to the excess mass spectrum from NA60 [21]. The actual 
NA60 acceptance filter and mass resolution have been incorporated [31]. The solid lines 
show the HSD results for a scenario including the collisional bradening of the p-meson 
whereas the dashed lines correspond to calculations with 'free' p spectral functions for 
reference. The dash-dotted lines represent the HSD calculations for the 'dropping mass + 
collisional broadening' model (see text). 

in Fig. [2] the corresponding NA60 data for the excess mass spectra with the HSD 
calculations for the dilepton yield from the direct p mesons, where the normalization 
to the data is performed with respect to the integral yield from 0.2 to 0.9 GeV for 
all scenarios considered. Note that an absolute normalization to the NA60 data - 
contrary to the systems discussed below - is not yet available and we thus focus on 
a 'shape' analysis. 

The dashed lines give the dilepton mass spectra for the direct p meson decays when 
incorporating only the vacuum p- spectral function in all hadronic reaction processes. 
These reference spectra overestimate the data in the region of the p-meson pole mass 
and underestimate the experimental spectra in the region below (and above) the pole 
mass such that in-medium modifications can clearly be identified. The solid lines 
show the result from HSD in the 'collisional broadening' scenario where no shift of 
the p pole mass is incorporated but an increase of the p-meson width due to hadronic 
collisions proportional to the baryon density ps (calculated as p 2 B {x) = j^(x)j^(x) 
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with denoting the baryon 4-current). The dash-dotted lines represent the HSD 
calculations for the 'dropping mass + collisional broadening' model where the p mass 
has been dropped with baryon density in accordance with Eq. (10) in Ref. |28j. An 
explicit representation of the p— and u— meson spectral functions employed here is 
presented in Fig. 2 of Ref. [28J. As can be seen from Fig. [2] the dilepton mass spec- 
trum is rather well described up to invariant masses of 0.9 GeV in the 'collisional 
broadening' scenario. For higher invariant masses the data signal additional contri- 
butions. This result is practically identical to the calculations of van Hees and Rapp 
[32] in the expanding fireball model - when incorporating the spectral function from 
Ref. p| - and demonstrates that the dominant in-medium effect seen in the p + p~ 
spectra from NA60 is a broadening of the p meson. The 'dropping mass + collisional 
broadening' model performs worse since it shifts too much strength to lower invariant 
dilepton masses. Note, however, that a reduced dropping of the vector-meson masses 
might be also compatible with the present data sets. 

In Ref. [33] a very detailed analysis of the NA60 data has been performed and shown 
that the additional yield above 0.9 GeV partly is due to open charm decays, four-pion 
collisions or 'quark- ant iquark' annihilation. We mention that our HSD calculations 
give only a small contribution from tt + a\ collisions in this invariant mass range but 
a preliminary study within the Parton-Hadron-String-Dynamics (PHSD) model [31] 
suggests that - apart from open charm decays - the extra yield seen experimentally by 
NA60 should be due to massive 'quark-antiquark' annihilations. Since this question 
is presently open and discussed controversely [3"5l36|37|l3"8"f 39] we concentrate on low 
mass dilepton pairs in the following. 

The next step in our study is related to an update of the HSD calculations in com- 
parison to the recent data from the CERES Collaboration [22] (with enhanced mass 
resolution). In Fig. [3] we present the HSD results for the mass differential dilepton 
spectra in 30% central Pb + Au collisions at 40 A GeV (upper part) and 7% cen- 
tral collisions at 158 A GeV (lower part) in comparison to the data from CERES 
[22] . The experimental dilepton yields in Fig. [3] are normalized to the average num- 
ber of charged particles for the corresponding centrality in the CERES acceptance: 
for 7% most central 158 A GeV events < N ch >= 177 and for 30% central 40 A 
GeV events < N c h >= 216 in the pseudo-rapidity interval 2.1 < n < 2.65. Cor- 
respondingly, the calculated dilepton yields have been also normalized to < N c h > 
obtained directly from the HSD calculations which are in good agreement with the 
measured numbers indicated above. The dashed lines in Fig. [3] show the results in 
case of vacuum spectral functions (for p,u,<f>) whereas the solid lines correspond to 
the 'collisional broadening' scenario. As in Fig. 2 the dash-dotted lines represent the 
HSD calculations for the 'dropping mass + collisional broadening' model. Similar to 
the In + In case the experimental data agree with the calculations employing the 
'collisional broadening' scenario (also in line with Ref. [33]). However, the CERES 
data also compare reasonably well with the 'dropping mass + collisional broadening' 
model up to invariant masses of 0.8 GeV. 
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Fig. 3. The HSD results for the mass differential dilepton spectra in 30% central Pb + Au 
collisions at 40 A- GeV (upper part) and 7% central collisions at 158 A GeV (lower part) 
in comparison to the data from CERES [22J. The dashed lines show the results for vacuum 
spectral functions (for p,u,(f>) whereas the solid lines correspond to the 'collisional broad- 
ening' scenario. The dash-dotted lines represent the HSD calculations for the 'dropping 
mass + collisional broadening' model (see text). 



On the other hand, the HSD model underpredicts the yield between the u and peaks 
which might again be attributed to possible contributions from 'quark- ant iquark' 
annihilations etc. We mention that in Ref. [18] we showed that this invariant mass 
region is very sensitive to the in-medium scenario since the simple 'dropping mass' 
picture provided a strong shift of the dilepton yield to the low mass regime and 
a strong reduction of the dilepton yield at M ~ 0.9 GeV whereas the 'collisional 
broadening' scenario gave a much larger contribution at M ~ 0.9 GeV. The 'drop- 
ping mass + collisional broadening' model (dot-dashed lines) shows a similar but 
less pronounced trend. Thus detailed measurements of this mass regime with high 
resolution will provide interesting constraints on the in-medium scenarios and might 
as well indicate possible contributions from partonic degrees of freedom - e.g. quark 
- antiquark annihilation or gluon-Compton scattering - already at SPS energies. A 
more detailed investigation within the PHSD approach will be presented in the near 
future. 
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Fig. 4. The comparison of the CERES data for the dilepton yield after subtraction of the 
hadronic cocktail (without the p) |22j with the HSD calculations: the dashed line shows 
the result for the vacuum p spectral function whereas the solid line corresponds to the 
'collisional broadening' scenario; the dash-dot-dot line stands for the excess yield defined 
as total yield after HSD 'cocktail' subtraction. The dash-dotted line represents the HSD 
calculations for the 'dropping mass + collisional broadening' model. 



In addition to the total yield shown in Fig.[3]the CERES Collaboration presented also 
the 'excess' yield (similar to NA60) denned by the dilepton yield after subtraction 
of the hadronic cocktail [22]. The hadronic 'cocktail' is composed by the sum of 
7!-°,r],uj,ri' Dalitz decays and direct decays of to and (f> mesons to dileptons. Thus, 
the residual dilepton yield might be attributed to the p meson contribution. The 
experimental excess yield is shown in Fig. H] (solid dots with errorbars) in comparison 
to the HSD calculations: the dashed line shows the result for the vacuum p spectral 
function, the solid line corresponds to the 'collisional broadening' scenario and the 
dash-dotted line represents the HSD calculations for the 'dropping mass + collisional 
broadening' model. The dash-dot-dot line stands for the excess yield defined by the 
total yield (for the 'collisional broadening' scenario) after subtraction of the HSD 
'cocktail'. As seen from Fig. HI in addition to the p meson contribution the HSD 
excess yield shows an enhancement at low invariant mass which is related to other 
dilepton sources not included in the 'cocktail', in particular to a sizeable contribution 
from the A Dalitz decay. The HSD excess yields show a reasonable agreement with 
the CERES data in the 'collisional broadening' scenario and in the 'dropping mass 
+ collisional broadening' model except of the first data point, which we relate more 
to uncertainties in the definition of the 'hadronic cocktail' than to new effects. More 
precise experimental data at very low mass as well as improvements of the 'cocktail' 
will shed more light on this issue. 

We step on to RHIC energies and first compare the HSD results for the dilepton 
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Fig. 5. The HSD results for the mass differential dilepton spectra in case of pp collisions 
at = 200 GeV m comparison to the data from PHENIX 00]. The actual PHENIX 
acceptance and mass resolution have been incorporated (see legend for the different color 
coding of the individual channels). 



invariant mass spectrum from pp collisions at y/s = 200 GeV with the data from 
PHENIX jlQ] in Fig. [51 The electron-positron pairs simulated in each pp event have 
been passed through the PHENIX acceptance and mass resolution routines [40J . We 
note that HSD for (inelastic) pp reactions for ^/s > 2.6 GeV is identical to FRITIOF 
(including PYTHIA v5.5 with JETSET v7.3 for the production and fragmentation 
of jets), i.e. for general hadron production. Explicit comparisons with data for pp 
reactions are presented in Refs. [H] as well as in Ref. [32] (Fig. 1). Since FRITIOF 
is modeled (fixed) to describe accurately elementary channels such as pp or tt + p 
reactions also HSD performs in a similar way. 

As seen from Fig. the HSD calculations well reproduce the PHENIX experimental 
spectrum which can entirely be described by meson Dalitz and direct decays as well 
as some contribution from open charm decays (light blue thin solid line as calculated 
by PYTHIA). This comparison demonstrates that the hadron production channels in 
HSD for elementary pp collisions are well under control also at the top RHIC energy. 

We recall that HSD also provides a reasonable description of hadron production in 
Au + Au collisions at ^fs = 200 GeV [JT] such that we can directly continue with the 
results for e + e~ pairs which are shown in Fig. [6]m case of inclusive Au + Au collisions 
in comparison to the data from PHENIX [23] . Again the actual PHENIX acceptance 
filter and mass resolution have been incorporated [43] . In the upper part of Fig. [6]the 
results are shown for vacuum spectral functions (for p, u, 0) including the channel 
decompositions (see legend for the different color coding of the individual channels). 
Whereas the total yield (upper solid line) is quite well described in the region of the 
pion Dalitz decay as well as the u and <fi mass regime we clearly underestimate the 
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Fig. 6. The HSD results for the mass differential dilepton spectra in case of inclusive 
Au + Au collisions at ^/s = 200 GeV in comparison to the data from PHENIX [23J. 
The actual PHENIX acceptance filter and mass resolution have been incorporated [33]. In 
the upper part the results are shown for vacuum spectral functions (for p,u,(f)) including 
the channel decompositions (see legend for the different color coding of the individual 
channels). The lower part shows a comparison for the total e + e~ mass spectrum in case of 
the 'free' scenario (dashed line), the 'collisional broadening' picture (solid line) as well as 
the 'dropping mass + collisional broadening' model (dash-dotted line). 

measured spectra in the regime from 0.2 to 0.6 GeV by an average factor of 3. 

When including the 'collisional broadening' scenario for the vector mesons we achieve 
the sum spectrum shown by the solid line in the lower part of Fig. [6] which is only 
slightly enhanced compared to the 'free' scenario (dashed line). Thus the question 
emerges if the PHENIX data might signal dropping vector meson masses? To answer 
this question we have performed also calculations in the 'dropping mass + collisional 
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broadening' model where the p and 00 masses have been dropped with baryon density 
in accordance with Eq. (10) in Ref. [28]. The respective HSD results are displayed in 
the lower part of Fig. E] by the dash-dotted line and indeed show a further enhance- 
ment of the dilepton yield which, however, is only small in the mass range 0.2 GeV 
< M < 0.4 GeV such that also this possibility has to be excluded in comparison to 
the PHENIX data. 

In summary we have used the off-shell version of the relativistic HSD transport model 
for the calculation of dilepton spectra from elementary as well as nucleus-nucleus col- 
lisions. Whereas the presently available dilepton data at SIS energies of 1 to 2 A-GeV 
(from the HADES Collaboration [H]) are well described in the 'collisional broaden- 
ing' scenario [28] this also holds for low mass dimuon data from In + In collisions at 
158 A- GeV (from the NA60 Collaboration) as well as the low mass dilepton spectra 
for Pb + Au collisions at 40 and 158 A -GeV (from the CERES Collaboration). The 
'dropping mass + collisional broadening' model for vector mesons seems compatible 
with the CERES data at SPS energies but does not perform well for the NA60 data. 
However, the low mass dilepton spectra from Au + Au collisions at RHIC (from 
the PHENIX Collaboration) are clearly underestimated in the invariant mass range 
from 0.2 to 0.6 GeV in the 'collisional broadening' scenario as well as in the 'drop- 
ping mass + collisional broadening' model, i.e. when assuming a shift of the vector 
meson mass poles with the baryon density. We mention that our results for the low 
mass dileptons are very close to the calculated spectra from van Hees and Rapp as 
well as Dusling and Zahed [45] (cf. the comparison in Ref. [46J). Consequently we 
attribute this additional low mass enhancement seen by PHENIX to non-hadronic 
sources, possibly to virtual gluon-Compton scattering [47] . 
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